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Abstract 

High quality jet noise spectral data measured at the Aeroacoustic Propulsion Laboratory at the NASA 
Glenn Research Center is used to examine a number of jet noise scaling laws. Configurations considered 
in the present study consist of convergent and convergent-divergent axisymmetric nozzles. Following the 
work of Viswanathan, velocity power factors are estimated using a least squares fit on spectral power 
density as a function of jet temperature and observer angle. The regression parameters are scrutinized for 
their uncertainty within the desired confidence margins. As an immediate application of the velocity 
power laws, spectral density in supersonic jets are decomposed into their respective components 
attributed to the jet mixing noise and broadband shock associated noise. Subsequent application of the 
least squares method on the shock power intensity shows that the latter also scales with some power of 

the shock parameter (M 2 . -M] ) . A modified shock parameter is defined in order to reduce the dependency 

of the regression factors on the nozzle design point within the uncertainty margins of the least squares 
method. 


Nomenclature 


A 

B 

c 

D 

f 

He= f D/Coo 
M=U/c 
M j ~U j 'cj 

M a =U/Coo 

m 

n 


P 


Pref 

PSD = lOLog[p 2 Uj/p^ ef D 
St= fD/Uj 


T 

U 


Ui 


Jet area 

Intercept parameter 

Sound speed 

Jet exit diameter (ft) 

Frequency (Hz) 

Helmholtz number 
Aerodynamic Mach number 
Ideally expanded Mach number 
Acoustic Mach number 
Shock intensity power factor 
Velocity power factor 
Acoustic pressure 

Acoustic reference pressure (0.0002 pbar) 

Power spectral density per Strouhal number, dB 
Strouhal number based on ideally expanded jet velocity 
Temperature ratio — normalized relative to ambient temperature 
Jet speed (fps) 

Ideally expanded jet speed (fps) 
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/ 3 Shock intensity parameter 

y Specific heat ratio 

6 Polar angle from jet inlet 

X Goodness factor 


Subscripts 

d 

e 

j 

s 

t 

oo 


Design point 
Nozzle exit 

Fully expanded condition 
Static condition 
Total or stagnation point 
Ambient conditions 


1.0 Introduction 

Development of a system level design and optimization tool is essential to the technology needs of 
the future commercial aircraft. As required by NASA’s Fundamental Aeronautics program, high fidelity 
design and analysis tools need to be developed to enable an optimal balance of performance, noise and 
emissions — and a rapid progression of new concepts and technology. Basic research is conducted within 
each discipline to gain an understating of the fundamentals, and to reduce uncertainty in modeling. 
Application of new technology to design efficient and environmentally friendly engines is an integral 
element of both Subsonic Fixed Wing and Supersonics projects. This article describes an ongoing activity 
in support of the environmental goals related to noise emission. 

Development of a successful analytical tool in jet noise prediction relies heavily on quality data. A 
typical narrowband spectrum in jet noise could consist of three distinct components, the jet mixing noise, 
shock-associated noise and screech. As the jet exit velocity becomes supersonic (relative to the ambient 
conditions), the mixing noise could be complemented by an additional component at small aft angles that 
is caused by the instabilities of the mean flow. From a modeling standpoint, it is extremely helpful if one 
succeeds in dividing the spectrum into components that are attributed to various generation mechanisms. 

The usual procedure is to start with a jet operating condition subject to a single generation 
mechanism. If one can propose scaling laws for this particular noise component, then it might be possible, 
with various degrees of success, to subtract that component from a multi-component spectrum. For 
example, shock-associated noise is obtained when we subtract the mixing noise from the total noise. Jet 
noise measurements are normally displayed as the power spectral density (PSD) per normalized 
frequency parameter, and are herein referred to as the total noise. When multiple sources are present, they 
are usually regarded as completely incoherent, and in the acoustic far field, the PSD of the sum is related 

, ,1 . 1 p.PSD / 10 r\PSD, no . 1 f\PSDj /10 

to the components as 10 =10 1 +10 2 

In aeroacoustics, the most quoted scaling law is the U s power law derived from Lighthill’s acoustic 
analogy (Ref. 1), as described in great detail in Reference 2. It concludes that the acoustic power emitted 
over the surface of a large sphere surrounding a jet is proportional to the eight-power of the jet exit 
velocity. Other attempts to fine-tune the power law (Refs. 3 to 5) consider the integrated power spectral 
density (i.e., the overall sound pressure level (OASPL) at each angle) to depend on both angle and jet 
temperature. Following the work of Viswanathan (Ref. 3), velocity power factors are evaluated from the 
narrowband data gathered at the Small Hot Jet Acoustic Rig (SHJAR) at the NASA Glenn Research 
Center, applying a least squares fit on the OASPL as a function of jet temperature and observer angle. The 
goodness of the fit is studied at each angle, and alternative relationships are proposed that improve the 
spectral collapse when certain conditions are met. 
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Since we use the integrated spectrum as the dependent variable in a linear regression, the PSD curve 
scales successfully from one jet to another only when (1) the basic shape of the spectrum is preserved, 
(2) the normalized frequency parameter meets the shift requirements. It will be noted that these two 
conditions are not always met at aft angles. However, it is possible to select a new normalized frequency 
parameter and collapse a major segment of the scaled spectra as long as the acoustic Mach number 
remains subsonic. 

Supersonic jets at imperfectly expanded condition emit broadband shock associated noise (BBSN) in 
addition to the jet mixing noise. In many cases, shock noise is accompanied with screech. Under some 
conditions (unheated jets in particular) the screech tones can be very strong and result in broadband 
amplification noise (AMN). This amounts to an additional 3 to 5 dB noise enhancement throughout the 
spectrum. Excitation of a jet by a pure tone has been established experimentally and there are indications 
that this amplification happens at sideline as well as small aft angles (Refs. 6 to 8). It is shown that 
scaling laws help identify the amplification. When intensity scaling is extended to BBSN, it makes a 
difference whether this amplification is added to the mixing noise prior to subtraction from the total noise, 
or if it is simply considered as an element of the shock noise. In the former approach, the mixing noise is 
increased uniformly and is matched with the total noise at the early part of the spectrum. In the latter 
approach, the unmatched mixing noise is subtracted from the total noise to obtain the BBSN. 

The paper is organized in the following order. Section 1 . 1 briefly describes the SHJAR experimental 
setup. In Section 2, we apply the linear least squares regression to find the velocity power factors for the 
jet mixing noise. The goodness of the fit, the joint confidence region of the two regression parameters, 
and their marginal inference interval is discussed at selective angles. Section 3 applies the scaling laws to 
a host of nozzles at various operating conditions and evaluates the quality of the spectral collapse at 
different angles and Mach numbers. The utility of the scaling laws are shown in Section 4 when we 
examine various noise components. In particular, the scaling laws are extended to the broadband shock 
associated noise and shock intensity parameters are examined versus angle and temperature. 

1.1 Experimental Setup 

Hot and cold jet noise spectral data have been acquired within the Small Hot Jet Acoustic Rig for an 
extensive test matrix. The SHJAR is a single stream hot jet rig located within the Aero Acoustic 
Propulsion Laboratory (AAPL) at the NASA Glenn Research Center. It uses remotely located 
compressors and a hydrogen burning combustor to provide heated air to test-nozzles at temperatures up to 
1300 °F and at pressures up to 150 psi. Additional noise generated by flow through piping and control 
valves is removed using a baffled muffler, and the flow is conditioned using screens to create clean and 
quiet conditions at the nozzle exit. Acoustic data repeatability has been determined to be within 0.50 dB 
in one-third octave bands at all frequencies (Refs. 9 to 12). 

The AAPL, which houses the SHJAR, is a 65-ft-geodesic dome lined with 24-in. -long sound 
absorbing wedges which remove sound reflection at all frequencies above 200-Hz and provide the 
anechoic environment required for noise studies. The floor and all surfaces around SHJAR are covered 
with fiberglass wedges (Fig. 1.1). The jet exhaust from SHJAR is directed outside through a large door. 
Additional details about the facility are provided in References 9 and 10. 

For most testing, SHJAR uses a 2.0 in nozzle, but can also operate larger nozzles with some 
limitation on cold set points at high Mach number. In typical testing, the SHJAR can cover a range of 
Mach numbers up to Mach 2.0, and static temperature ratios up to 2.8 using a hydrogen combustor and 
central air compressor facilities. Stagnation flow conditions for the rig are measured within a 24 in. 
plenum upstream of the nozzle. 

Nozzles used in the present noise measurements are designated as smcxxx, where “xxx” is a 3 digit 
number that identifies a particular nozzle (Table 1.1). The smcOOO is a baseline convergent nozzle with a 
5° conic contraction and a 0.04-in.-thick lip (Fig. 1.2), and smc021 is a modified smcOOO geometry in 
order to minimize the screech noise (Fig. 1.3). In addition to the data repeatability, SHJAR measurements 
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Figure 1.1. — Small Hot Jet Acoustic Rig (SHJAR). 



TABLE 1.1 NOZZLE CONFIGURATIONS TESTED 
IN SHJAR FACILITY 


Nozzle 

Configuration 

M d 

D e - in. 

NPR d 

smcOOO 

Convergent 

1.00 

2.0 

1.89 

smc021 a 

Convergent 

1.00 

2.0 

1.89 

smc014 

CD 

1.185 

2.0 

2.37 

smc015 

CD 

1.40 

2.0 

3.18 

smc016 

CD 

1.50 

2.0 

3.67 

smc018 

CD 

1.80 

2.0 

5.74 


Convergent nozzle with reduced screech. 



Figure 1.3. — smcOOO and smc021 nozzles. 


have also been examined versus other existing measurements (Refs. 3, 13 to 17) in order to identify 
possible extraneous sources that are unrelated to jet noise. 

Measurements are based on a test matrix designed to isolate various parameters that are thought to 
drive the physics of noise generation and its propagation such as constant jet velocity and constant 
temperature (static or stagnation). For jets operating at a pressure ratio larger than 1.89, i.e., ideal Mach 
number >1.0, convergent-divergent nozzles were designed for ideal expansion and were tested at different 
temperatures. The result is a matrix of data at fixed acoustic Mach number M a , aerodynamic Mach 
number M and temperature ratio T for a broad range of conditions, with and without screech and 
broadband shock-associated noise. 


2.0 Mixing Noise and Velocity Power Factors 

Power laws for the jet mixing noise are derived from a matrix of test cases that are free of shock- 
associated noise. A convergent nozzle (i.e., smcOOO, Table 1.1) operating at a nozzle pressure ratio (NPR) 
<1.89 is considered for this purpose (NPR is defined as the total plenum pressure divided by the ambient 
static pressure). A set of N test cases is tabulated at a constant static or stagnation temperature (Ref. 5), 
where each case corresponds to an exit velocity £/., i = 1,2,..., A. Measurements are reported on a 100- 
diameter (200 in) arc, and at 24 equally spaced angles 6 (50° to 165° relative to the nozzle inlet), and 
have been corrected for atmospheric attenuation. 
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The integrated spectrum (i.e., the overall sound pressure level OASPL) at each jet velocity U i is 

y t = OASPL{6J\ i = 1,2 ,...,7V (2.1) 


The integration is carried out in a narrowband, and within a fixed Strouhal frequency range of 
0.02 < St < 10.0 after removing possible extraneous facility noise at either low or high end of the 
spectrum. Using a linear least squares method, we define a line-fit 

y, = n(0,T)x. +B(0,T); x = 10 Log^ / cj, i = 1,2, ...,7V (2.2) 

n and B are referred to as the velocity power factor and the intercept parameter respectively. The intercept 
parameter B is the same as the approximated overall sound pressure level directivity (OASPL) at 
temperature T if the exit acoustic Mach number were 1.0. Subsequently the effect of the arc distance on 
noise is included in parameter B. As pointed out earlier, all measurements presented here are at an arc 
distance of 100D and are corrected for the atmospheric attenuation. 

Let the uncertainty at each data point be <T , then the measure of goodness of the fit is defined as 


mr) = 


1 f iy-y ,) 2 

N-2tt a, 


(2.3) 


We require TV to be at least 3. Additionally, all experimental data are considered to have equal 
uncertainty (<7 = 1) . Other statistical properties of the linear fit such as confidence intervals for 
regression parameters n and B and their joint confidence region (Ref. 18) are also examined at selective 
angles. 

When the regression errors are within acceptable bounds, we extrapolate the scaling laws to the jet 
mixing noise from supersonic jets. A Scaled PSD is defined once we remove the extra dB associated U" 


from each jet 


scaled PSD = PSD - \0n(6,T)Log(U j /cJ-10Log(A. / A e ) (2.4) 


The scaled spectra are now expected to collapse in amplitude across different jets, including 
supersonic jets. As stated earlier, the normalized frequency parameter and the basic shape of the PSD 
curve need to meet certain requirements for the collapse to be successful. The last term on the right hand 
side of Eq. 2.4 signifies that in imperfectly expanded supersonic jets, a fully expanded velocity U. is 

reached at an area A. that is different from the exit area A e . Overexpanded jets reach full expansion 

prior to the nozzle exit, while underexpand jets reach a fully expanded jet velocity downstream of the exit 
plane. When supersonic jets are operating at off-design conditions, the mismatch between the exit static 
pressure p e and the ambient pressure p ^ results in the shock formation. If the associated entropy change 

is considered small enough to permit an isentropic flow relation pip 7 - const , then the fully expanded 
area A is estimated as 
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p, 

p. 


(2.5) 


f l + M 7 2 (f-l)/2V-‘ 
ll + M 2 (r-l)/2j 

y + 1 

A. M,( p\*r 
A e M.\p J 

where M . is the fully expanded Mach number, and M d is the design Mach number. In addition, the usual one- 
dimensional compressible flow relations hold between the nozzle pressure ratio NPR, Mach number M . and the 
corresponding static temperature T s 


y — 1 — 

npr = a + - — M 2 y-' . 
y 2 J> 

r-i 

TIT = NPR 7 . 

t s 


( 2 . 6 ) 

(2.7) 


Using the convergent nozzle smcOOO, a host of SHJAR test points at subsonic conditions (i.e., NPR 
<1.89) were grouped according to static temperature T s = 1.0, 1.20, 1.43, 1.76, 2.27, 2.70. Within the first 
two groups (i.e., T s = 1.0, 1.20), jets were confined to subsonic acoustic Mach numbers, while the latter 
four groups consisted of jets at both subsonic and supersonic acoustic Mach numbers. The result of the 
least squares fit is usually displayed as a set of parametric curves for regression parameters n(0, T ) and 
B(0,T ) as in Figure 2.1. The outcome is expected to be sensitive to the number of test points in each 
group and the range of the dependent variable X- . Detailed listings of the set points and the sensitivity 
analysis of the regression parameters are provided in Reference 5. Figure 2.2(a) shows the goodness 



Figure 2.1. — Velocity power factor n and the intercept parameter B as a function of angle and static temperature. 
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Figure 2.2 — (a) Goodness factor for the regression parameters of Figure 2.1 . (b) Goodness of the fit if the regression 
were limited to points at subsonic acoustic Mach number only. 



Figure 2.3 — (a) Measurements (symbols); least squares fitted line (red); and the 95% confidence band (dashed 
lines), (b) Joint 95% confidence region (ellipsoid) and marginal inference intervals (dashed lines) of parameters 
nand B.(T S = 2.70,0 = 60°) 


factor corresponding to Figure 2.1. A small goodness factor is indicative of a successful regression. When 

data variance is relatively large (as in Fig 2.2(a) ^ - 120 ), the addition of new set points to the group 
does not necessarily improve the linear regression. However, if set points were to be confined to M a < 1, 
then data scatter would remain small at all angles as shown in Figure 2.2(b). 

To examine the confidence band of the interpolation parameters, consider the measurements at 
T = 2.70, where 7 points (0.39 < M< 1.47 ) have been utilized in the regression. Figure 2.3(a) shows the 
regression points, the least squares fit (red line) and the 95 percent confidence margin at 60 . The joint 
confidence region (Ref. 18) falls within n = 4.89 ±0.17 and B = 92.56 ± 0.33 as seen in Figure 2.3(b). The 
error band is usually expected to decrease with increasing the number of regression points. However, this 
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Figure 2.4— (a) Measurements (symbols); least squares fitted line (red); and the 95% confidence band (dashed 
lines), (b) Joint 95% confidence region (ellipsoid) and marginal inference intervals (dashed lines) of parameters 
nand B. (T s = 2.70, 0= 160°) 


is not true at small aft angles where the range of the dependent variable x. plays a more instrumental role 

due to the increased scatter at M a >1. For example at <9=150°, six interpolation points with M a <1 
generates a much better fit compared to 8 points at both M a < 1 and M a >1 . 

Shown in Figure 2.4 is the 95 percent confidence band of the two parameters at 160°. The joint 
confidence region is now increased to n- 6.87 ± 0.70, and B = 107.82 ± 1.48 . The above error estimates 
point to the deteriorating quality of the linear regression method at aft angles. It will be shown that 
changes in the spectral shape with jet speed at aft angles make it increasingly difficult to obtain a 
satisfactory spectral collapse throughout the frequency range at all speeds. However, the regression 
parameters n and B in Figure 2.1 could still be employed to collapse, at least, a major segment of the 
subsonic data at aft angles. 

3.0 Examination of Jet Noise Spectral Data Using Velocity Power Factors 

In order to equip us with the necessary tools to examine the jet noise power spectral laws, additional 
charts for the regression parameters n and B are evaluated as a function of angle and stagnation 
temperature. Figure 3.1 shows such a template when interpolation consisted of 5 points with subsonic 
aerodynamic Mach number at each temperature. Details of the error analysis with and without points at 
supersonic acoustic Mach number are given in Reference 5. Conclusions are similar to those stated in the 
previous section and show a gradual decline in the quality of the regression fit at aft angles. Since 
parameter B measures the OASPL directivity at M a = 1.0 where the two temperatures convert as 
T t = T s + ( y - 1) /2 , then a one-to-one relationship exists between the two intercept parameters at related 
temperatures. A parallel correspondence is not possible between the two n factors because velocity cannot 
be held constant in order to convert temperatures. 

According to Figure 3.1(b), parameter B appears to collapse at the three higher stagnation 
temperatures of 1.8, 2.2 and 2.7 — subsequently velocity scaling of jet noise across such temperatures is 
achieved with a single parameter. Lower temperature jets, on the other hand, show slight decrease in 
parameter B with increasing temperature (Fig. 2.1(b) at static temperatures of 1.0 and 1.20), subsequently 
both parameters n and B would be required for a successful scaling regardless of which temperature is 
used. The major difference between the two temperatures relates to approximations made in defining a 
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Inlet angle, 0 


Inlet angle, 0 


Figure 3.1 . — Velocity power factor n and the intercept parameter B as a function of angle and stagnation temperature. 


single static temperature when in reality it is affected by changes in the velocity profile. Experimentally, 
the plenum total temperature is a control parameter, while the static temperature at the jet exit is 
calculated using one-dimensional Bernoulli’s equation for a compressible flow. 

We now inspect the spectral collapse at selected angles, temperatures, and Mach numbers using a 
combination of different nozzles geometries. 

The operation conditions for 9 jets at a constant stagnation temperature ratio of 1.0 are shown in 
Table 3.1. The first seven readings are due to a convergent nozzle smcOOO (see Table 1.1 for nozzle 
geometry), with two readings 1616 and 1618 at underexpanded conditions. Reading 1636 is due to a CD 
nozzle operating at its design condition, i.e., M d =1. 50, NPR d =3.671 . The last reading (1605) is a 
convergent nozzle smc021, however the nozzle lip is modified with fine notches in order to minimize the 
screech-related noise at supersonic conditions. 


TABLE 3.1.— SHJAR READINGS AT PLENUM TEMPERATURE RATIO 1.0 


Rdg 

Nozzle 

T 

1 s 

T, 

U./c 

J 00 

M 

NPR 

M. 

J 

A. /A 

J e 

1610 

smcOOO 

0.97 

1.0 

0.40 

0.40 

1.117 

0.40 

1.0 

1611 


0.96 

1.0 

0.49 

0.50 

1.186 

0.50 

1.0 

1612 


0.91 

1.0 

0.67 

0.70 

1.387 

0.70 

1.0 

1613 


0.86 

1.0 

0.83 

0.90 

1.692 

0.90 

1.0 

1614 


0.83 

1.0 

0.91 

1.00 

1.893 

1.00 

1.0 

1616 


0.76 

1.0 

1.08 

1.00 

2.556 

1.24 

1.043 

1618 


0.70 

1.0 

1.23 

1.00 

3.514 

1.47 

1.156 

1636 

smc016 

0.69 

1.0 

1.24 

1.50 

3.671 

1.50 

1.0 

1605 

smc02 1 

0.80 

1.0 

1.04 

1.00 

2.328 

1.17 

1.02 


Figure 3.2 shows the 90° scaled power spectral density for all seven convergent nozzles following 
Equation 2.4. With the exception of the underexpanded cases, all spectra collapse at n = 7.93 (see Fig. 3.1 
for power factors) and within the uncertainty of the regression parameters. The two underexpanded jets at 
readings 1616 and 1618 contain shock-associated noise as well as intense screech. The screech tones 
appear as spikes in the spectra, and if one succeeds in eliminating such tones, then the jet mixing-noise is 
expected to dominate the spectra at low frequency even in the presence of the shock-associated noise. 
This is illustrated in Figure 3.2(b) with convergent nozzle smc021 at reading 1605, as well as 
Mach 1.50 CD nozzle at reading 1636. The lack of collapse of underexpanded spectral data to the level of 
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Figure 3.2. — Application of velocity power law at conditions of Table 3.1 (7) = 1 .0, 0 = 90°, n = 7.93). 
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Figure 3.3. — Application of velocity power law at conditions of Table 3.1 . (a) 7) = 1 .0, 0 = 50°, n = 8.1 1 . (b) 7) = 1 .0, 
0= 120°, n = 8.30. 
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a fully expanded jet at low frequency, as seen in jet reading 1618, is attributed to an amplification of 
broadband jet noise due to screech. An amplification of -4.0 dB is noted in this case. 

Experimental measurements of the turbulent kinetic energy in jets with and without screech show an 
increase in turbulence level with screech. Recent PIV measurements of Bridges et al. (Ref. 19) indicate a 
direct correspondence between broadband amplification noise and the additional turbulence measured in 
the presence of screech. Shear-layer excitation at a correct frequency often explains the reason for 
differences in noise measurements between different rigs. Jet noise prediction models that use a RANS 
input are unable to capture this amplification due to the inability of a time-averaged calculation to account 
for tone excitations. Such amplifications need to be removed from data prior to validation of predictions. 

The scaled PSD is additionally examined at 50 and 120° in Figure 3.3 using the appropriate velocity 
power factor n(0,T) . As the observer progressively moves towards aft angles, the BBSN noise weakens 
gradually, however, amplification effects may remain equally strong at all angles. 
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3.1 Collapse of Spectral Data at Small Aft Angles 

Two fundamental assumptions are buried in the spectral power law (Eq. 2.4). In the first place, it is 
implied that the basic shape of the spectrum at any angle 6 , described by the first term on the right, is 
invariant with respect to jet velocity and temperature. Second, since the spectral density is expressed as a 
function of the jet Strouhal frequency, spectra are assumed to peak at the same Strouhal number. When 
any of these assumptions are violated, data will collapse poorly, and other measures are necessary to fix 
that. 

Figure 3.4(a) illustrates an attempt to collapse the subsonic noise data at 150° (see Table 3.1). An 
inspection of the spectral data confirms that neither of the rules cited above hold at this angle. First, the 
spectra appear relatively wider at lower jet speeds, and display a broader peak. Higher speed jets, on the 
other hand, exhibit a well-defined peak. It is readily demonstrated that removing the lowest speed jet from 
the group would increase the value of the power factor n from 9.13 to 9.32. Removing the second lowest 
speed jet will still increase n to 9.72. This is expected, as the area under the spectrum (on a Fog basis) is 
highly dominated by the peak. When the peak covers a wider span, that particular jet would have an 
exaggerated effect on the power factor. The goodness factor described in the previous section remains 
small only when the general shape of the spectrum is preserved. Second, the peak Strouhal number is seen 
to vary between 0.18 and 0.35. 

Both problems are fixed once we increment the velocity power factor by 1.0, i.e., from 9.13 to 10.13, 
and select the Helmholtz number (fD/cJ) as the normalized frequency parameter as seen in 
Figure 3.4(b). 

To see this more clearly, note that the power spectral density per Helmholtz number converts to that 
per Strouhal number according to 


10 Log(p 2 Coo / p 2 ref D ) = 10 Log(p 2 U. / p 2 ref D ) - 10 Log(U / cj (3.1) 


Therefore, the left hand side of Equation 2.4 may also be interpreted as the scaled power spectral 
density per Helmholtz number provided that n is replaced with (n+\) on the right hand side. For brevity, 
we choose “scaled PSD ” to represent either of the two meanings depending on the selection of the 
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Figure 3.4. — Application of the power law at conditions of Table 3.1 (7)= 1.0, 0 = 150°). (a) Scaled power spectral 
density per Strouhal no. (n = 9.13). (b) Scaled PSD per Helmholtz no. (n = 10.13). 
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Figure 3.5. — Application of power law. (a) smcOOO and smc016 nozzles at set points of Table 3.1 (7) = 1 .0, 0 = 150°, 
n = 9.13+1). (b) smcOOO and Mach 1.18 smc014-CD nozzles at indicated Mach numbers (7) = 1.80, 0 = 160°, 
a? = 8.90+1). 


frequency parameter (and subsequent adjustment in parameter n). The difference in spectral shapes is now 
clearly seen as an increase in the low frequency noise in lower jet speeds, while the spectral amplitudes at 
higher frequencies collapse together for all subsonic conditions. 

Failure of the velocity scaling of supersonic jet noise spectra at 150° is seen in Figure 3.5(a) when a 
shock- free CD nozzle (reading 1636 — Table 3. 1) is added to the subsonic jets. A better presentation of the 
supersonic spectra at small aft angles is shown in Figure 3.5(b) at 160° and at heated condition of 
T t = 1.80. The first 5 jets are due the convergent nozzle smcOOO, and the last jet is smc014-CD nozzle at 
its design point M = 1.18. The scaled spectra show that the low frequency noise approaches a limit with 
increasing jet speed. As the acoustic Mach number exceeds 1.0, a broadening takes place at the spectral 
peak, and the high frequency roll off fails to collapse with a single velocity power factor. The high 
frequency noise shows a collapse as long as the acoustic Mach number remains subsonic. If the collapsed 
segment of the subsonic spectra is attributed to a particular generation mechanism, then it could arguably 
be scaled up to the velocity of a supersonic jet and be viewed as one component of noise at such speeds. 
The supersonic jet mixing noise at aft angles may, by itself, be considered as superposition of two 
incoherent components. More detailed study of the spectral features of jet noise at small aft angles is 
provided in Reference 5. 


4.0 Shock Noise and Intensity Power Laws 

As pointed out previously, a typical jet noise spectrum could be a composite of three distinct 
components, i.e., jet mixing noise, shock-associated noise, and screech. Shock noise is due to the interaction 
of the jet turbulence with the shock cell structures under imperfectly expanded conditions, and dominates 
the spectra at mid to high frequency, mostly at forward and sideline angles (Refs. 20 to 23). Harper-Bourne 
and Fisher (Ref. 24) pioneered a formalism that estimates the shock noise component of supersonic jets. 
Screech noise, on the other hand, appears as isolated spikes at almost all angles when a jet is imperfectly 
expanded. Strong jet-related screech tones often excite the instabilities of the jet and enhance the broadband 
jet noise. The so-called amplification noise (AMN) may be identified when we scale up the jet mixing noise 
from a shock- free subsonic condition to the speed of a supersonic jet at the same temperature 

PSD 2 =PSD l +\Qn(P 9 T)LogQJ j2 I U jX ) + \0Log(A j2 / A jX ) (4.1) 
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Subscripts 1 and 2 refer to the lower-speed (i.e., datum) and higher-speed (i.e., destination) jets 
respectively. The lower speed jet is required to be shock-free (i.e., A jX !A e = 1). Let the spectral density 

associated with the mixing noise component in the higher speed jet be PSD 2 . This component is evaluated 
by up scaling shock-free subsonic spectra at a similar angle and temperature as per velocity power laws. 
When amplification noise is minimal, PSD 2 should coincide with the total spectral density of jet2 at the 
initial segment of the spectrum (usually at St < 0.10), otherwise the difference is attributed to the 
amplification effect. To isolate the shock-associated noise in the higher speed jet, we subtract the mixing 
noise from the total noise. When amplification is not negligible, it may be accounted for two ways: 

(1) added to the mixing noise to calculate the matched-mixing noise or, (2) it could be included in the 
shock noise by subtracting the unmatched mixing noise from the total. Both methods are examined here. 

A systematic approach to isolating the shock noise is presented in the following order: 

(1) An under expanded convergent nozzle without broadband amplification effect AMN 

(2) Under expanded convergent nozzle with intense screech and amplification AMN 

(3) Convergent-divergent (CD) nozzles. 

Within the first category a reduced-screech convergent nozzle at a nozzle pressure ratio >1.89 would 
be appropriate. An example is configuration smc021-1605 from Table 3.1 at NPR = 2.32 (fully expanded 
Mach number of 1.17). To isolate the mixing noise component, velocity scaling is applied using any one 
of the subsonic jets at a similar stagnation temperature T t = 1.0 in Table 3.1. Figure 4.1(a) shows the total 
noise (blue), and the mixing noise (dark) at inlet angles of 50 to 120°. A nearly perfect match at the low- 
frequency part of the spectra points to the absence of the amplification effect. Broadband shock associated 
noise BBSN is subsequently evaluated as 

^ qPSD (shock)/ 10 y qPSD (total)/ 10 ^qPSD (mixing)/ 10 (4 2 ^) 



Figure 4.1. — Underexpanded convergent nozzle smc021 at reading 1605 (Table 3.1). (a) Total noise (blue); jet 
mixing noise (dark), (b) Shock-associated noise. 
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Shown in Figure 4.1(b) is the BBSN, with peak values of 1 12.5dB at (St = 0.90, 6 = 50°), and 103. 3dB at (St = 2.2, 
6 = 120 ). The primary peak frequency of the BBSN is given as 


/,= 


C 

L{\ + M c cos#) 


(4.3) 


where L- 1.10-^Af 2 - 1 D according to HBF formulation (Ref. 24), and L-(n / o x )^M y 2 - 1 D. following 


Tam’s formulations (Refs. 20 and 25) with cr l =2.40 , and M c is the convection Mach number of the 
turbulent eddies U c /c co . 

The shock-noise intensity, i.e., integrated spectrum in Figure 4.1(b), usually scales as J3 m , where 
shock parameter jS= ^ M 2 . - M ] is a measure of deviation from the design point, and m is the shock 


power factor. A nominal value for m in chocked convergent nozzles is 4.0. Experimental measurements 
(Ref. 26) show that m varies with both angle and temperature and approaches 4.0 in the very forward 
angle where BBSN is at its peak. Here we use the SHJAR measurements to examine a similar intensity 
law for the shock-associated noise. The linear regression method of Section 2 is applied as: 


y i =OASPL shock (0,n i = 1,2,..., A 

y t = m(0,T)x i + B(0,T); x f = 10 Log ft , 



(4.4) 


The intercept parameter B is a measure of the BBSN intensity (dB) when J3 = 1. In convergent nozzles 
(Md = 1.0) this occurs at a nozzle pressure ratio of y rl(r ~ l) . As before, all intensity calculations are 
reported at 100D arc, and as lossless. In order to examine the intensity scaling for BBSN, two nozzles are 
considered here — a reduced screech convergent nozzle smc021 at underexpanded conditions of Table 4.1, 
and the smc016 CD nozzle at both over- and under expanded conditions of Table 4.2. 


TABLE 4.1— TEST CONDITIONS— CONVERGENT 
NOZZLE SMC021 GROUPED AT CONSTANT 
VALUES OF STAGNATION TEMPERATURE 


Rdg 

T, 

U./c 

J 00 

NPR 

M. 

J 

A. /A 

J e 

2079 

1.0 

0.95 

2.00 

1.048 

1.00 

2080 

1.0 

1.04 

2.36 

1.178 

1.02 

2081 

1.0 

1.18 

3.15 

1.392 

1.11 

2082 

1.0 

1.24 

3.64 

1.493 

1.17 

2083 

1.0 

1.34 

4.65 

1.659 

1.30 

2085 

1.0 

1.40 

5.70 

1.795 

1.43 

2088 

1.80 

1.26 

2.00 

1.049 

1.00 

2089 

1.80 

1.38 

2.35 

1.181 

1.02 

2090 

1.80 

1.57 

3.15 

1.391 

1.11 

2091 

1.80 

1.65 

3.63 

1.492 

1.17 

2092 

1.80 

1.77 

4.64 

1.659 

1.30 

2093 

1.80 

1.85 

5.69 

1.791 

1.43 

2100 

2.20 

1.39 

2.00 

1.045 

1.00 

2098 

2.20 

1.53 

2.35 

1.180 

1.02 

2097 

2.20 

1.73 

3.14 

1.390 

1.11 

2096 

2.20 

1.82 

3.62 

1.493 

1.17 

2095 

2.20 

1.96 

4.64 

1.659 

1.30 

2094 

2.20 

2.05 

5.67 

1.793 

1.43 

2102 

2.70 

1.54 

1.99 

1.048 

1.00 

2103 

2.70 

1.69 

2.33 

1.174 

1.02 

2104 

2.70 

1.92 

3.12 

1.390 

1.11 

2105 

2.70 

2.02 

3.60 

1.489 

1.17 

2106 

2.70 

2.17 

4.63 

1.660 

1.30 

2107 

2.70 

2.27 

5.65 

1.789 

1.43 
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TABLE 4.2— TEST CONDITIONS— CD NOZZLE 
SMC016 (M d = 1.50) AT STAGNATION 
TEMPERATURE RATIO T t = 1.0 


Rdg 

U./c 

J 00 

NPR 

M. 

J 

A. /A 

J e 

1636 

1.246 

3.669 

1.500 

1.0 

1637 

1.250 

3.701 

1.505 

1.002 

1638 

1.254 

3.739 

1.512 

1.006 

1639 

1.266 

3.851 

1.532 

1.016 

1640 

1.285 

4.046 

1.566 

1.033 

1641 

1.309 

4.323 

1.611 

1.058 

1642 

1.339 

4.704 

1.667 

1.092 

1643 

1.372 

5.207 

1.735 

1.137 

1644 

1.418 

5.995 

1.827 

1.206 

1628 

0.947 

2.003 

1.047 

0.889 

1629 

1.051 

2.405 

1.193 

0.901 

1630 

1.119 

2.749 

1.294 

0.923 

1631 

1.174 

3.098 

1.381 

0.950 

1632 

1.212 

3.383 

1.443 

0.974 

1633 

1.226 

3.498 

1.466 

0.984 

1634 

1.238 

3.602 

1.486 

0.993 

1635 

1.244 

3.654 

1.496 

0.998 




Figure 4.2. — Decomposition of total noise into mixing and shock noise components in Mach 1.49 convergent nozzle 
smc021 at NPR = 3.64 (reading 2082, Table 4.1). Total noise (blue); jet mixing noise (dark); shock-associated 
noise (red), (a) Mixing noise excludes amplification, (b) Mixing noise includes amplification. 

The convergent nozzle is examined first. Figure 4.2(a) illustrates the application of the velocity power 
law to a typical jet condition (i.e., reading 2082 — Table 4.1) whereby the total noise is decomposed into 
its two components. Unlike the previous example shown in Figure 4.1, the mixing noise and the total 
noise fail to collapse together at the early stages of the spectra. The difference (-2.0 to 3.0 dB) is due to 
the amplification noise (AMN). Shown in Figure 4.2(b) are the matched-mixing noise and the subsequent 
BBSN. Note that the shock-associated noise is slightly different depending on how the amplification is 
book kept. 


NAS A/TM— 2009-2 15674 


15 


6 




Figure 4.3 — (a) BBSN versus shock parameter (3 in the convergent nozzle smc021 at Tf = 1 .0, 0 = 50°: data (sym- 
bols); least squares fitted line (red); and the 80% confidence band (dashed lines), (b) 80% Joint confidence margin 
(shaded region) and marginal inference intervals (dashed lines) for shock power factor m and intercept B. 




Inlet angle, 0 Inlet angle, 0 

Figure 4.4. — Shock factors versus angle and temperature in an underexpanded convergent nozzle smc021 when 
amplification is included in mixing noise, (a) Shock power factor m(d,Tf). (b) Intercept parameter B(Q,Tf). 

Application of the least squares method to the BBSN at <9= 50° and T t = 1.0 at moderately 
underexpanded conditions is shown in Figure 4.3. Set point 2085 at NPR = 5.70 has been removed from 
least squared fit in order to reduce the data scatter and improve the confidence margins. In this figure, the 
AMN noise component has been added to the mixing noise similar to Figure 4.2(b). At 80 percent 
confidence margin, the interpolation parameters cover inference intervals of m- 3.89 ± 0.90 , 
B =115.93 ± 2.15 . 

Figure 4.4 shows the two parameters m(0,T t ) and B(0,T) for the convergent nozzle at four 
temperatures. The unheated jets show a distinctly different set of parameters, and at a constant M radiate 
less BBSN at forward angles where shock noise its relatively strong. Figure 4.5 shows the shock 
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Inlet angle, 0 Inlet angle, 0 

Figure 4.5. — Shock factors versus angle and temperature in an underexpanded convergent nozzle smc021 when 
amplification is included in shock noise, (a) Shock power factor m(0,7"^). (b) Intercept parameter B(0,7))- 




Figure 4.6. — Noise components in convergent nozzle smc021 and CD nozzle smc016 versus fully expanded 
Mach number Mj at 50° and at 7) = 1 .0. (a) Amplification is added to mixing noise, (b) Amplification is included 
in BBSN noise. 


parameters for the same nozzle conditions if amplification noise were considered as part of the BBSN and 
component decomposition were carried out according to Figure 4.2(a). As before, the heated jets display 
distinctly different characteristics with their amplitude parameter B clustered together for most of the 
angle range. The marginal inference intervals of parameters m and B (not shown here) are quite 
comparable irrespective of the decomposition method. 

Next, we examine the shock parameters in the Mach 1.50 CD nozzle (smc016). Test conditions 
consist of unheated set points only, as listed in Table 4.2, including both under and over expanded 
conditions. As before, when screech-related amplification is present, noise elements are distinguished 

using two methods of component splitting. Sample result at <9 =50° is shown in Figure 4.6, where 
amplification is added to the mixing noise (part a), or is considered as part of the shock noise (Part b). The 
same figure also displays noise components in the convergent nozzle smc021 following a similar 
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Figure 4.7 — BBSN versus shock parameter p in the CD nozzle smc016 at = 1 .50 and 7) = 1 .0, 0 = 50°: data 
(symbols); least squares fitted line (red); 95% confidence band (dashed lines), (a) Underexpanded regime. 

(b) Overexpanded regime. 




Inlet angle, 0 Inlet angle, 0 

Figure 4.8. — Shock noise power factors A77(0,T)) and intercept factor B(0,7)) in the CD nozzle smc016. 

(a) Amplification is added to mixing noise, (b) Amplification is included in the BBSN noise. 

decomposition. Point P3 displays the design conditions for the CD nozzle. Away from the design point, 
BBSN intensifies in both under expanded (P3-P4) and over expanded (P3-P2) regimes of the CD nozzle. 

A linear interpolation of the BBSN is relatively successful only when we confine the least squares 
method to the moderately imperfectly expanded conditions, subsequently interval points (P2-P1) are 
excluded from the following linear interpolation. In highly overexpanded jets, a normal shock occurs 
within the nozzle and flow decelerates through the remaining diverging section and exhausts as a 
subsonic jet. Experimental results of (Ref. 27) indicate that an additional broadband noise component 
could also be generated in the low Mach number range of the overexpanded regime, depending on the 
expansion angle of the divergent section. 

Figure 4.7 shows sample interpolations at 0=50° with readings (1637-1644) and (1630-1634) from 
Table 4.2, and following component splitting of Figure 4.6(a). The joint confidence margins are 
m = 2.17 ± 0.32, B = 121.17 ± 1.50 in the underexpanded regime, and m = 2.45 ± 0.77, B = 116.31 ± 2.49 in 

the overexpanded regimes, respectively. 

Shown in Figure 4.8 are the BBSN intensity interpolation parameters m(0,T ) and B(6, T ) in the 
unheated CD nozzle using two different methods for bookkeeping of the amplification noise. The shock 
power factor m is slightly different in the under and over expanded regimes, however within the 
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uncertainty of the regression m ~ 2.0 in the forward angles, which is almost half of that for a convergent 
nozzle at a similar temperature (Figs. 4.4 and 4.5). 

4.1 Modified Shock Parameter 

Interpolation parameters m(6,T) and B{6, T ) associated with the BBSN as evaluated in the previous 
section depend on the nozzle design point. For example, in the convergent nozzle with unheated jets, 
m = 4 and B = 1 17.0 dB at forward angles (Figs. 4.4 and 4.5) — while in the CD nozzle m = 2.0 and B = 
116.0 dB or 121. OdB depending on the overexpanded or underexpanded regimes, respectively (Fig. 4.8). 
It is practically more desirable to obtain a set of parameters that are independent of the nozzle design 
point or mode of operation in the CD nozzles. Subsequently, a modified shock intensity parameter /3 m is 
defined as 


y 6 = (j8 t M 3 J 2 ) d , underexpanded 

(4.5) 

= (A, j‘' » overexpanded 

Only CD nozzles are affected by the new definition. As seen in Figure 4.9, there is a relative collapse 
of the intercept parameter B between different modes of operation. In addition, parameter B is also in 
agreement with the convergent nozzle at a similar temperature, and within the uncertainty of regression 
scheme. The new estimates for the shock power factor m also exhibit less discrepancy relative to the 
convergent nozzle at aft angles. The gap in power factors at aft angles remains as before due to a steady 
reduction in m with increasing angle# in overexpanded jets. Also shown in Figure 4.10 are typical 
regression results using the new parameter (3 m . The rescaled shock parameter affects the horizontal 
distribution of regression points as shown in Figures 4.10(a) and 4.7(a). 
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Figure 4.9. — Shock noise power factors m(d,Tt) an d intercept factor B(0,Td in the CD nozzle smc016 with a modified 
shock parameter as defined in Equation 4.5. (a) Amplification is added to mixing noise, (b) Amplification is 
included in the BBSN noise. 
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Figure 4.10 — (a) BBSN versus modified shock parameter (3 m in the convergent nozzle smc016 at 7) = 1.0, 0 = 50°: 
data (symbols); least squares fitted line (red); and the 95% confidence band (dashed lines), (b) 95% Joint confidence 
margin (shaded region) and marginal inference intervals (dashed lines) for shock power factor m and intercept B. 
(Amplification is added to the mixing noise). 


5. Summary 

An extensive database of jet noise has been analyzed, correlating power spectral density as a function 
of angle and temperature over a range of jet velocities in single stream jets. Velocity scaling of jet mixing 
noise was shown to be equally successful at constant static or stagnation temperatures. In general, two 
parameters are needed to achieve a successful scaling when jets are at different temperatures, although at 
higher temperatures, a single parameter n will suffice at most angles as shown in Figure 3.1(b). Spectral 
collapse of jet noise data at small aft angles is less successful due to changes in the spectral shape with jet 
speed. A collapse of the high frequency roll off of amplitudes at aft angles is achieved with Helmholtz 
frequency only when the acoustic Mach number is subsonic. 

Velocity scaling of the jet mixing noise helps identify other noise components, such as broadband 
amplification of jet noise in the presence of screech, and the broadband shock associated noise. Two 
different methods were discussed in bookkeeping of the amplification noise, and subsequent shock noise 
intensity parameters were evaluated for convergent and convergent-divergent nozzles using a linear 
regression method. Finally, a modified shock intensity parameter was introduced that allows the two 
regression parameters to be less dependent on the nozzle design point within the uncertainty margins of 
the interpolation method. 

Scaling laws should help with jet noise data evaluation, component identification, and most 
importantly with the development of jet noise prediction tools. 
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